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Abstract. We investigate the relationship between persistent currents in multi-channel rings containing an 
embedded scatterer and the conductance through the same scatterer attached to leads. The case of two 
uncoupled channels corresponds to a Hubbard chain, for which the one-dimensional embedding method 
is readily generalized. Various tests are carried out to validate this new procedure, and the conductance 
of short one-dimensional Hubbard chains attached to perfect leads is computed for different system sizes 
and interaction strengths. In the case of two coupled channels the conductance can be obtained from a 
statistical analysis of the persistent current or by reducing the multi-channel scattering problem to several 
single-channel setups. 
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1 Introduction 

Electronic correlations can influence coherent electronic 
transport in striking ways. An example is the Kondo effect 
in the transport through quantum dots [1-3]. It is well- 
known that the inclusion of the many-body effects arising 
from electron-electron interactions in quantum transport 
calculations is extremely difficult. When a coherent and in- 
teracting nanosystem is attached to Fermi liquid reservoirs 
via leads, the fundamental problem consists in the match- 
ing of the correlated many-body wave-function in the 
nanosystem with the effective one-body wave-functions in 
the leads that allow to define the bias voltage as difference 
between the chemical potentials of the reservoirs. Even 
for the simplest nanosystems such a matching procedure 
is very complicated to apply [4-6]. There is an increasing 
literature where the auxiliary Kohn-Sham wave- functions 
of density functional calculations are used to determine 
the transport properties of realistic systems. However, no 
rigorous theoretical justification for this approach exists 
at present and approximations in the density functional 
can considerably influence the results for transport [7]. 

At sufficiently low temperatures and when the inter- 
actions are restricted to a small region in space that is 



attached to noninteracting leads, a Landauer-like formula 
with an effective interaction-dependent transmission [8] 
holds for the linear (dimcnsionlcss) conductance, 



.9 = Tr{tt^}. 



(1) 



For leads supporting N open channels, the transmission 
matrix t is of dimension N. Determining this transmission 
matrix remains however a hard task. 

The recently developed embedding method [9-15] rep- 
resents a considerable advance for the single-channel case, 
N = 1, allowing to extract the effective transmission prob- 
ability \t\ 2 of the nanosystem at the Fermi energy of the 
reservoirs. This method is based on the fact that the per- 
sistent current through a ring, made of a scatterer and a 
non-interacting lead and threaded by a magnetic flux <f>, 
depends on the transmission probability of the scatterer. 
At zero temperature the flux dependence of the ground- 
state energy E determines the persistent current through 



dE 

' d<j> 



(2) 



In the limit of very large; lead length Ll and for a non- 
interacting scatterer, the persistent current and the trans- 
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mission probability can be easily related. For an odd num- 
ber of spinless fermions in the ring, the flux dependence of 
the persistent current in lowest order in l/L as a function 
of the transmission amplitude t at the Fermi wavenumber 
kp is given by [16] 



Ttx\ ev F arccos(|i(fc F )|cos(<2>)) . 
yfi - \t(k F )\ 2 cos 2 (@) 



(3) 



Here, L$ denotes the length of the nanosystem, L = L$ + 
Ll is the total length of the ring, and <P = (e/h)<fi the 
dimcnsionless flux. 

A similar expression holds for an even number of fer- 
mions. In particular, 



Iff 



l-oo V J° L (n/2) 



(4) 



where J®(ir/2) is the persistent current of a perfect ring 
without scatterer at flux <P — tt/2. 

Using the same relationships for the interacting case 
allows to extract \t\ from the many-body ground state 
properties of the ring. It has been shown [14] that for 
an interacting single-channel scattering problem the full 
flux dependence of the persistent current can be described 
by d3l). Such a check provides a strong evidence for the 
validity of the embedding method. A nanosytem can then 
be replaced by an effective one-body scatterer, provided 
the electrodes are connected to the interacting region by 
sufficiently long noninteracting leads [15, 17]. 

A convenient measure for the persistent current is the 
phase sensitivity V = (L/2)\E ^ — 2?o,a|) where E ^ P and 
.Eo,A denote the ground state energies of the ring for peri- 
odic and anti-periodic boundary conditions, respectively. 
An extrapolation of T> for L — > oo keeping Lg constant 
yields 2?oo and the absolute value of the effective trans- 
mission amplitude of the system at the Fermi energy 



t = sin 



7T £>oc 

2VL 



(•5) 



Here, = ttHvf/4, with the Fermi velocity wp, is the 
phase sensitivity of a ring with perfect transmission t = 1 
through the scatterer. 

Another quantity that can be used to characterize the 
flux dependence of the persistent current is the curvature 



C = L 



e2 d^E 



(6) 



<p=o 



of the many-body ground state energy as a function of 
the flux. Because of a level crossing occurring at zero flux, 
the curvature diverges for the case of an even number of 
particles in a clean ring. For an odd number of particles, 
the curvature is well-behaved. Its long- lead limit L — *■ oo 



r — 

*-T>0 



e 2 VF \t\ arccos(|i|) 



(7) 



can be obtained from j3J, and depends only on the abso- 
lute value of the transmission amplitude [18]. The inverse 



of ([7]) is unique, but in general \t\ has to be determined 
numerically from . 

It is important to stress that the previous relations be- 
tween the persistent current and the transmission proba- 
bility at the Fermi energy have only been established for 
the case of spinless fermions with strictly one-dimensional 
leads. According to (fTl, the conductance in the case of 
more than one channel is given as the sum over all trans- 
mission probabilities between pairs of channels. Then, the 
value of the persistent current alone is in general no longer 
sufficient to determine the conductance. 

In the present paper we generalize the embedding 
method to fermions with spin and multi-channel situa- 
tions. After a discussion of persistent currents in rings 
embedding a multi-channel scatterer (Sec.pl, we address 
the one-dimensional Hubbard chain attached to perfect 
leads (Sec. [3]). In this latter case, because of the spin rota- 
tion symmetry it is straightforward to extend the embed- 
ding method developed for the single-channel case in order 
to include spin. We discuss some technical aspects of the 
evaluation of the conductance, develop an improved em- 
bedding method based on damped boundary conditions, 
and present numerical results for short chains. 

Two paths are presented in order to treat the case of 
two coupled channels. The first one is a statistical study 
of the persistent currents with random channel mixing 
(Sec. [i]). The second one is a reduction of the multi- 
channel problem to several single-channel problems, al- 
lowing to infer the conductance from the corresponding 
single-channel transmission coefficients that can be de- 
termined usingthe standard procedure of the embedding 
method (Sec. [5] and Appendix [C) . This second path is 
used in order to argue that channel mixing is absent in 
the case of Hubbard chains with spin-rotation symmetry. 
We have relegated to the appendices some of the lengthy 
formulations and alternative proofs. 



2 Persistent currents in multi-channel rings 
embedding a scatterer 

In this section we address the relationship between the 
transmissions of a multi-channel non-interacting scatterer 
and the persistent current in rings embedding the scat- 
terer. In particular, we underline the difficulties in going 
from one to two channels. 



2.1 AT-channel scattering problem 

The ideal quasi-one-dimensional leads are characterized by 
translational invariance along the wire axis and a position- 
independent lateral confinement, allowing to separate the 
one-body Hamiltonian into a longitudinal and a trans- 
verse part. The solutions of the longitudinal part are plane 
waves and the eigenstates 4> n (y) of the transverse part 
give rise to the conduction channels and energy offsets e n . 
It is convenient to describe the states in the lead at a 
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Fig. 1. Sketch of the transfer-matrix description of scattering 
in multi-channel wires. 



given energy e as a superposition of flux-normalized one- 
dimensional plane waves as 



iV 



&(x, y) = [A n exp (ik n x) + B n exp (-ik n x)] <f) n (y) 

n=l v fc " 

(8) 

where k n denotes the wave vector in the nth channel. In 
general, the channel index n also accounts for the elec- 
tron spin as in the case of the Hubbard chain discussed in 
Section [3] 

As represented in Figure [TJ the transfer matrix Ms 
relates the amplitudes A n and B n of W\ on the left-hand 
side to the corresponding amplitudes C n and D n of W\i on 
the right-hand side of the scatterer according to 



(9) 



Assuming current conservation and time reversal sym- 
metry, it is possible to express the transfer matrix of an 
A^-channel scatterer in the polar decomposition 







MA 


Cn 




A N 


D x 


= M s 


Bi 


\D N ) 




\B N ) 



urv 



(10) 



as a product of three matrices, where 



U = 
V = 

r = 



u 

u* 

v 

v* 



(11) 



VT+x Vx 
Vx ViTx 



u and v are complex unitary N x iV-matrices, and A is a 
real, diagonal N x N- matrix [19]. The matrix A contains 
all transmission probabilities |i a | 2 = 1/(1 + A a ) for the 
different eigenmodes of the scatterer, where A a is the a th 
entry of A. Therefore, A determines the conductance of the 
scatterer through g — J2a=i V(l + Aq ) . The matrices u 
and v describe the way these eigenmodes are connected to 
the different incoming and outgoing channels. In general, 
a mixing of channels as illustrated in Figure [2] will occur. 
The main interest of the polar decomposition consists in 





u 


r 
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Fig. 2. A general transfer matrix Ms satisfying current con- 
servation and time-reversal symmetry can be expressed as the 
product ( |10[ ) of three matrices. U and V relate the eigenmodes 
in the leads to those in the system and f describes the trans- 
mission and reflection properties of the system. 



the very different role played by the radial (A) and an- 
gular parameters (u,v), allowing the development of ran- 
dom matrix theories to describe the quantum transport 
through an arbitrary scatterer [20,21]. 

An iV-channel ring consisting of a scatterer described 
by Ms and an ideal lead can be characterized by the to- 
tal transfer matrix of the ring Ms Ml, where Ml is the 
diagonal transfer matrix of the lead with entries 



(Ml)„„ 

( M L)«+jV,n+AT 



exp(r# + ifc„£ L ), 
exp (i<£ - ifc„i L ) , 



(12) 



describing right- and left-moving electrons (1 < n < N). 

The one-particle eigenstates ip( pS> of the ring are de- 
fined by the amplitudes A^ and through the matrix 



condition 



0. 



(1-M S M L ) 

The corresponding eigenenergies e^ p ' satisfy 
det(l - M S M L ) =0, 



(13) 



(14) 



where 1 denotes the 2A^ x 2A^ identity matrix. The con- 
tribution j W to the persistent current of the ring follows 
from the flux dependence of through Equation (14). 



Alternatively, it can also be obtained from the expectation 
value of the current operator in the leads in the limit of 
very large L L , 



i(p) _ 



eh 



N 



\B 



(P)|2 



n=l 



(15) 



Summing the contributions of all occupied states p we 
are lead to the persistent current through the ring. 



2.2 One-channel scattering problem 

For energies e between the two lowest transverse energies 
E\ and £2 only one channel propagates and the problem 
is particularly simple [14] since u and v are pure phases 
given, respectively, by exp[i(a±C)/2] in terms of the angles 
a and £. The quantization condition ( 14 1 can be written 
as 

cos(^) = cos(fcL + 8a) . (16) 



sin(#) sin(y) 



4 



A. Freyn et al: The embedding method beyond the single-channel case 



Following the standard parametrization of a 2 x 2 
transfer matrix [14], we have defined the angles 9 and ip 
through 

'1 + Acos 2 (0 



sin(0) = 
sin(yj) = 



1 + A 
1 



(17) 



Vl + Acos 2 (C) 



In this way t = exp(ia) sin(0) sin(y) is the transmission 
amplitude and Sa — a — kLg is the scattering phase- shift. 
We remark in passing that the expressions simplify in the 
particular case with right-left symmetry, where ip = ( = 
tt/2. 

T akin g the flux derivative for the p-th solution of Equa- 
tion ( 16 I and following the line of reasoning presented in 



the appendix of [14] we obtain the contribution which 
after summation over the occupied states leads to pj). For 
the particular case of <P = tt/2 we have [14] 



j (p) (7r/2) = (-l) p ~ fc (p) |t(/c (p) )| . 

IjTTh 



(18) 



The alternative route through the eigenvectors is also 
simple in the N = 1 case since Equation (15) becomes 



with 

lx (p) l 2 



B 



dp) 



(v) 



eh 
mL 



i-|x (p) l 2 



Ap)\ 2 



(19) 



(p) 



1 



1 - sin 2 (0) sin 2 (<^) 
x [l + sin 2 (0)sin 2 (<^) - 2 sin 2 (0) sm 2 {ip) cos 2 (<£)- 
2 sin(0) sin((/?) sin(#) 
x sin [±arccos(sin(0) sin(<p) cos(^))] ] , 

(20) 

where the ± refers to the parity of the eigenstate label p. 
Using 

>)|2 = (lTsin(6>)sin(y)) 2 
1 - sin 2 (0) sin 2 (v?) 



IX 1 



l<f=ir/2 



(21) 



where X t = rf/t 2 , with < ti,n < 1 and t 2 +r| = 1. The 
mixing angle ip introduced here should not be confused 
with the angle ip introduced above for the one-channel 
case. 



The quantization condition ( 14 1 can now be written as 
2cos 2 (0) -Tr{M}cos(<£) + F(M) -1=0, (23) 



where M = M s M L (0 = 0) and F(M) is given by the 2 x 2 
sub-determinants of M as 



F(M) — 

3=2 



Mh My 

Mn Mjj 



(24) 



Assuming a parabolic dispersion relation in the leads, 
ki and k-i are related by h 2 (k\ — k 2 )/2m = e 2 — £i — Ae. 
Defining k = (fci + fc 2 )/2 we have k\ — fc 2 = mAe/h 2 k 
and the quantization condition (23) at <P = 7r/2 can be 
expressed as 



F(k, L L , t u t 2 , if, a, /3, 7, ip, e, r), ( 



(25) 



The allowed k^ define the eigenstates of the ring, with 
associated persistent currents 



(tt/2) 



ehk^ /Tr{M}\ 
to [dF/dk] 



(26) 



fc=fc(p) 



In Appendix |A| we give the expressions of Tr{M} and 
^(M) in terms of k\ ± k%, Ll, and the angles char- 
acterizing the matrices u and v. The complexity of such 
expressions forces us to introduce several simplifications 
in our problem. The most radical one is to take the two 
channels as spin-degenerate modes. This approximation 
is treated in the next section, where we analyze the Hub- 
bard chain and show that most of the complexities of the 
two-channel problem are not present in this case. In Sec- 
tion [4] we consider a simplified two-channel problem and 
introduce statistical concepts to extract the transmission 
coefficients. 



in Equation (pi we are lead to the contribution (pi to 3 Conductance of a Hubbard chain 



the persistent current. 



2.3 Two-channel scattering problem 



For energies e satisfying e 2 < e < £3, we have two propa- 
gating channels and the matrices u, v and A appearing in 
(11) can be parametrized as 



v = 



A = 



cos(<p)e ia sin(v3)e i ( Q+ ' 3 )\ 
- sin(^)e i7 cosMe^^ ) ' 

i cos(-0)e 1£ — isin(^)e 
-i sin(i/;)e i ( £+ '') -i cos(V>)e i(e+, ' : 

Ai 
A 2 



(22) 



In this section, we consider electron transport in a chain 
with Hubbard- like electron-electron interactions in a finite 
segment of the chain. This is a particular case of a two- 
channel chain where the spin degree of freedom of the 
electrons gives rise to the two channels. We will show that 
the embedding method can be applied in this special case 
in a quite straightforward manner and precise numerical 
calculations can be carried out. 



3.1 Model, conductance and spin-rotation symmetry 

The Hamiltonian of the whole system reads 

H = H K + H V , (27) 
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where 



E E( c I,^ +1 , CT +h.c.) (28) 



i <7=T,1 



is the homogeneous kinetic energy part describing elec- 
trons in an ideal ID chain. Here, Ci, CT annihilates an elec- 
tron with spin a on site i, and we define the energy scale 
by setting the hopping amplitude equal to one. The in- 
teracting region situated on sites 1 to i§ is distinguished 
from the rest of the chain solely by the presence of the 
on-site Hubbard interaction 



ffu = C/E(^T-l/2)(n U -l/2) 



(29) 



with strength U. Here, hi <a = c\ cr c i(J . Thus, our inter- 
acting system is a chain of L$ sites on which a spin-up 
electron interacts with a spin-down electron on the same 
site. 

We fix the Fermi energy in the leads to the center 
of the band, corresponding to half-filling. The interaction 
term H\j contains a background potential which renders 
the system particle-hole symmetric and ensures half-filling 
in the interacting region independent of the interaction 
strength. 

The two channels in our problem are defined by the 
two possible orientations of the electron spin a — {t, !}• 
In the absence of a magnetic field, the corresponding states 
in the leads are degenerate. 

If we are allowed to describe the many-body scatter- 
ing through the interacting segment of the chain by an 
effective single-electron scattering situation, the conduc- 
tance is given by the sum over all effective transmission 
probabilities as 

9 = ^ Ka'f , 



(30) 



CT,<T> 



where f CT)Cr / is the transmission amplitude from the incom- 
ing channel of spin a into the outgoing channel character- 
ized by spin a' . 

We show in the sequel that spin-rotation symmetry im- 
poses serious limitations on the transmission coefficients. 
Since the Hamiltonian is invariant under a rotation of the 
spin basis {|T), II)}, there can be no preferential spin ori- 
entation. As a consequence, the transmission matrix must 
have the same symmetry. Invariance of the transmission 
matrix t under arbitrary spin rotations excludes transmis- 
sions accompanied by a spin flip so that 



0. 



(31) 



Otherwise the transmission matrix could be diagonalized 
by a spin rotation thereby violating the spin-rotation sym- 
metry. Furthermore, the transmissions with spin conser- 
vation have to be equal 



IT 



(32) 



An alternative way of demonstrating the conditions 



(31) and (32) is to use the procedure of the reduction to 
the single-cnannel scattering. We develop this second path 
in Appendix |B| 



3.2 Numerical results for the conductance of Hubbard 
chains 

As discussed in the previous section, the effective one- 
body scattering describing the transmission through a 
Hubbard chain is characterized by spin conservation and 
the two spin channels are not mixed by the effective one- 
body scatterer. The transmission through the effective 
scatterer is therefore characterized by a single parame- 
ter, the spin-independent transmission amplitude t. The 
conductance ( 30 1 then simplifies to 



ff=2|t| 



(33) 



This expression holds for arbitrary interaction strength 
in our Hubbard chain. The spin appears as a factor of 
two, as in the case of general one-dimensional Fermi liquid 
systems with spin [13]. 

The relation between persistent current and effective 
one-body transmission amplitude which forms the basis 
of the standard embedding method for the single-channel 
case can now be generalized to two channels without mix- 
ing. The persistent current is then given by the sum of 
the contributions from the two independent channels. If 
in addition to the spin-rotation symmetry of the Hamilto- 
nian, we choose a spin- independent filling of the two sub- 
systems with iVj = JVj_, the two contributions are equal 
such that the two-channel persistent current is twice the 
single-channel value. As a consequence, §5§ is replaced by 



\t\ = sin 



4 V" 



(34) 



(2) 

where T>ca is the numerically obtained two-channel phase 
sensitivity in the limit of a very large ring. When the ring 
size is not much larger than the Hubbard chain, the flux 
dependence of the current may change considerably and 
reach a very different functional form for Ll <C [22] . 

We use the Density Matrix Renormalization Group al- 
gorithm (DMRG) [23-25] adapted to the Hubbard model 
to calculate T>^> for a two-channel ring with a Hubbard 
segment, taking fully into account the electronic correla- 
tions. These many-body ground state properties are cal- 
culated for half-filled rings with increasing sizes L corre- 
sponding to odd numbers — = L/2 of electrons 
for each spin channel. We keep up to 6000 states in the 
DMRG iterations, and perform fits to extrapolate the size- 
dependence of £>( 2 ) to [18]. 

The numerical results for the conductance of short 
Hubbard chains using (34 1 are presented in Figure [3] For 
even L$ we obtain a suppression of the conductance by the 
interaction that becomes more pronounced for longer Hub- 
bard chains. This can be interpreted as a precursor of the 
Mott-Hubbard transition. For weak interaction strengths, 
our numerical results presented in Figure [3] are in agree- 
ment with the results of the perturbative approach of 
Reference [26]. However, the perturbative approach does 
not yield the significant independence that we obtain 
at strong interaction. Nevertheless, a strong interaction- 
induced reduction of the conductance has been obtained 
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Fig. 3. Conductance of short Hubbard chains as a function of 
the interaction strength. The solid lines are guide to the eyes. 
The dashed curve shows the perturbative result produced using 
the theory of Reference [26]. 



for a Hubbard chain of length Lg with reduced coupling 
to the reservoirs [27]. 

For odd values of Lg it has been predicted from pertur- 
bative arguments [26,28] and a numerical renormalization 
group (NRG) study [29] that the transmission should be 
perfect. As shown in Figure [3] this perfect conductance 
g = 2 is found in the case of a single Hubbard impurity 
Lg = 1. For other odd numbers of interacting Hubbard 
sites, our results at strong interaction show small devia- 
tions from g — 2 (see Table [I] for the Lg = 3 case) . These 
deviations are of the order of the uncertainty of the ex- 
trapolation procedure included in the embedding method. 
It is then important to develop complementary approaches 
to test, and eventually improve, these values. This task is 
undertaken in the next two subsections. 



3.3 Numerical limitations of the standard embedding 
method 

As mentioned in the introduction, the curvature yields an 
alternative route for extracting the absolute value of the 
transmission amplitude. In the two-channel case, Equa- 
tion ^ becomes 



2e 2 vp \t\ arccos(|f|) 

00 = 



(35) 



In order to verify that the many-body two-channel 
scatterer modelled by a chain with Hubbard-like electron- 
electron interactions can indeed be replaced by the effec- 
tive one-body scatterer without channel mixing described 
in the previous section, we compute the transmission from 

the curvature Coo [30] and compare it with the one ob- 

(2) 

tained from the phase sensitivity L>oo . Since the latter is a 
measure for the integral of the persistent current over the 
flux interval < # < w, and the former corresponds to the 
slope of the increase of the persistent current as a function 

(2) (2) 

of flux at <P = 0, 2?oo and Coo characterize very different 
aspects of the persistent current and its dependence on the 



Table 1. Values of the effective transmission probability |t| 2 
for a Hubbard chain of length Ls = 3 obtained from the em- 
bedding method using the phase sensitivity (first line) and the 
curvature (second line), for different values of the interaction 
strength [7. The last line is the result of a two-parameter fit to 
the flux dependence obtained using damped boundary condi- 
tions. 





[7 = 




[7 = 4 


[7 = 6 


from T) 


1.00 


1.00 


0.98 


0.75 


from C 


1.00 


0.99 


0.95 


0.65 


damped boundary 




0.9996 


0.996 


0.968 



2.5 



1.5 



0.5 




tt/2 tt 

Fig. 4. Persistent currents extrapolated to infinite lead length 
as a function of the magnetic flux for a ring containing a Hub- 
bard chain with Ls = 3 and [7 = 6 (squares). The error bars 
indicate the difference between a linear and a parabolic extrap- 
olation of the available data. The line is a fit of (J3j) to the data 



points yielding \t\ 



0.87. 



flux. The two quantities are related to the transmission of 
an effective one-body scatterer by the relations (34 1 and 



(351. If the many-body scatterer can be characterized by 
an effective one-body scatterer, then the flux dependence 
of the persistent current in the correlated system has to 
agree with the one of a ring with an effective one-body 
scatterer, and the two alternative ways of extracting \t\ 
should be equivalent. This is the case for even Lg, where 
the results obtained from (35) are very close to the ones 
resulting from the phase sensitivity shown in Figure [3] 
The strongest deviations are of the order of 2 %, which 
is smaller than the symbol size. For odd Lg, the results 
of Table [T] show a difference between the two methods for 
large values of the interaction strength. However, these 
differences are of the order of the deviation from the ex- 
pected result of perfect transmission. Therefore, we can 
conclude that the two methods agree within their preci- 
sion. The consistency of the values of |t| 2 extracted from 

^ and can then be considered as strong evidence 
for the validity of the embedding method. 

For the most critical case of Table [l] (U = 6) we have 
performed the numerically very demanding calculation of 
the persistent currents for different values of the flux. The 
numerical data extrapolated to infinite lead length are 
shown in Figure [4] The uncertainty of the extrapolations 
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indicated by the errorbars remains unfortunately rather 
important. The one-parameter fit of ^ (with an addi- 
tional factor of 2 for the spin) to the data points yields 
\t\ 2 w 0.87, much closer to the expected perfect transmis- 
sion than the values extracted from T>^ and C^} shown 
in Table Q] 

In comparison with the single-channel case, the cor- 
related two-channel scattering problem requires consider- 
ably higher numerical efforts at equal ring length. As a 
consequence, we had to limit our computations to total 
ring sizes L < 30. As discussed above, this is not generally 
a problem, but it affects the quality of the extrapolations 
to infinite ring size for situations close to conductance res- 
onances [14], such as the case of half filling and odd values 
of Lg. The reliability of the extrapolations of the numeri- 
cally obtained values to infinite lead length becomes lim- 
ited when the energy resolution given by the level spac- 
ing in the largest rings is insufficient to resolve the reso- 
nance at the band center, whose width strongly decreases 
with increasing odd Ls and U. In order to overcome this 
problem, considerably larger ring sizes would have to be 
considered to improve the extrapolations. However, it ap- 
pears impossible within the standard embedding method, 
and with the present computational resources, to perform 
precise enough calculations for sufficiently large ring sizes. 
In the next subsection we develop an improved embedding 
method that allows to overcome these difficulties. 



3.4 Improved embedding method: Damped boundary 
conditions 

In order to overcome the limitation in ring size, we intro- 
duce damped boundary conditions [31] that allow to sim- 
ulate longer effective lead lengths. The damped boundary 
conditions consist in a reduction of the last Lrj hopping 
elements of the Hamiltonian on both sides of the Hubbard 
chain, leading to a region with reduced local bandwidth 
in the ring opposite to the Hubbard part, as indicated in 
Figure [5] by the distance between the upper and the lower 
line. The increased local density of states is reminiscent 
of a long lead having full hopping values. We choose an 
exponential reduction controlled by a damping parameter 
A < 1 such that the last hopping elements on either side 
of the interacting region are given by A 3 with j counting 
the last I/d sites. 

We have checked that for the cases of Figure [3] 
the results using damped boundary conditions reproduce 
those from the reliably extrapolated standard embedding 
method. However, the usefulness of the improved method 
appears in the cases previously discussed where the stan- 
dard embedding method meets its numerical limitations. 
We therefore show in Figure [6] numerical results for the 
persistent currents using damped boundary conditions for 
L§ = 10, with L-o = 10 and A — 0.7. The lines are fits 
of the analytical result for the persistent current of nonin- 
teracting particles without damped boundary conditions 
(formula (3| with an additional factor of 2 for the spin). 
The transmission of the scatterer and the effective length 



dam P ed boundary bandwidth 

-L D - 



interacting 
region 




Fig. 5. Sketch of the damped boundary condition setup. The 
hopping matrix elements are exponentially reduced on the Ld 
last sites on both sides of the interacting region, locally sup- 
pressing the bandwidth as indicated by the distance between 
the upper and lower line. The increased density of states in the 
damped boundary region simulates an effectively longer lead. 
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Fig. 6. Persistent currents as a function of flux for half-filled 
Hubbard rings with L — 30 and Ls = 3. Data points are 
obtained using damped boundary conditions with Ld = 10 
and damping A = 0.7. Crosses, pluses, and squares are for 
interaction strengths U = 2, 4, and 6, respectively. The lines 
are two-parameter fits of the formula Q yielding the values of 
Table [l] for the transmission and effective lengths L c g — 278, 
282, and 295, respectively. 



of the ring have been used as fitting parameters. The qual- 
ity of the fits is a strong indication for the validity of the 
embedding method in the Hubbard model case since it 
provides a comparison of the flux dependence with that 
of the noninteracting case as it was done for the single- 
channel case [14]. 

The fitting parameters show indeed that the effective 
lead length can increase considerably with the damping, 
provided it is not chosen too strong. Damping parame- 
ters larger than A — 0.7 result, for the system param- 
eters chosen here, in backscattering, and therefore ren- 
der the reliability of the extracted conductances question- 
able. Remarkably, the resulting transmission approaches 
the expected perfect conductance (see Table [l]) . Thereby 
we have obtained a considerable improvement from the 
standard extrapolation procedure, in particular for cases 
where the latter yields poor results. The results of the im- 
proved embedding method can thus be considered to be 
consistent with the perfect conductance for all odd values 
of is that is expected from perturbative arguments [26,28] 
and NRG [29]. 
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4 Persistent currents for two degenerate 
channels 

In the previous section we considered the case of de- 
generate and uncoupled channels, where the generaliza- 
tion of the one-dimensional embedding method is rather 
straightforward. Once we take into account channel mix- 
ing, the one-to-one correspondence between persistent cur- 
rents and transmission coefficients is lost. However, we still 
expect that the A parameters of the polar decomposition 
( 10 ) will be determining the corresponding persistent cur- 
rents. We will treat in this section an arbitrary channel 
mixing, and in order to keep the problem tractable we 
consider quasi-degenerate channels, where \ki — k 2 \ k. 
In this case, if we neglect the weak fc-dependence of u, v, 
and A, the solutions of (25 1 are doubly degenerate, since 



the roles of and A. 
changeable. 



(or and B 2 P> ) are inter- 



4.1 Symmetric nanosystems with fixed channel mixing 

In the case where the nanosystem exhibits left-right in- 
version s ymm etry we can set a — e, (3 — r\, 7 = 9, and 
(p = ip in ( 22 1 . This assumption considerably simplifies the 
expression of Tr{M} and _F(M) presented in Appendix |A| 
and thus the resulting analytical calculations. 

Even in the symmetric case the problem is still rather 
complicated. We therefo re co nsider particular values of ip. 
For instance, for ip — 0, (25 1 leads to 



cos(2[fcL L + a + P + 7]) = cos(2[a - /? - 7]) (36) 
and one gets two families of solutions 



k {e ' 1) L L =-2a + £TT 
k^L L =-2(f3 + 1 ) + en. 



(37) 



Using ( 26 1 , these quantized k- values lead to the single- 



levels currents 



eh 



^.») = (_!)< k W\t b \ 6 = 1,2, (38) 

Choosing (p = tt/2, one gets similar results with an in- 
terchange of ti and t 2 , and a relative sign between the 
two current contributions. The fact that the persistent 
currents associated with the ring eigenmodes are solely 
determined by the transmission eigenvalues is to be ex- 
pected since for these particular values of ip there is no 
mode- mixing by M. 

A more interesting case is that of tp = 7r/4, which corre- 
sponds to the maximum mode-mixing and the most likely 
value if we take the matrix u as uniformly distributed in 
the unitary ensemble [20] . The corresponding quantization 
condition is 



(-D b 

-a- /? - 7 + ^— ^arccos(x) + £tt (39) 



with 



X = Z\sin 2 (a - 7) + cos 2 (a - 7) cos(2/3) , 
A = txt 2 - r x r 2 . 



(40) 
(41) 



The resulting persistent currents are given by 

) 

--(-lY +l " Zr sgn[cos(a - 7)] 



3 



,b) _( ■ 

2mL^ 



(* 2 - ti)sgn[cos(/3)] 



l+A 



tan 2 (a — 7) 



+(-1) 



2cos 2 (/3) 

b (h + £ 2 )sgn[sin(/?)] 



(42) 



1 



l-A 
2sin 2 (/3) 



tan 2 (a — 7) 



and therefore are bounded by 



;(<,&) 



< — max{ii,i 2 j- 



(43) 



In the limit of large we can take k^ :1 ^ ~ fc^' 2 ) ~ fcW 
£tt/Ll and show that 



3^-3™ 



< — = — l*i — *a| , 



< 



(*i + *a) 



(44) 
(45) 



We remark that, leaving aside the trivial case of un- 
coupled channels, the relationship between persistent cur- 
rent and transmission amplitudes is not unique, as it in- 
volves the mixing angles. However, already the special case 
ip = 7r/4 shows that there are stringent constraints on the 
possible values of t x and t 2 . If we disposed of an ensemble 
of different scatterers with fixed ii, t 2 , and mixing an- 
gle ip, we could extract the transmission amplitudes by 
varying the angles a, (3 and 7 through the ensemble. We 
carry on this procedure by choosing these three angles uni- 
formly distributed in the interval (0,2-7r), solving the re- 
sulting eigenvalue problem (13 1, and then obtaining from 
( 15 1 the two resulting persistent currents. For ip — 7r/4 we 
obtain the distribution of normalized persistent currents 



J 



(1,2) 



= mLi^/(eh\k\)j^ 1 ' 2S> shown in Figure [7] where we 
have arbitrarily taken t\ — 0.8 and t 2 = 0.4. The restric- 
tions imposed by (44 1 and ( 45 ) constrain the possible val- 
ues of jW) to the rectangle defined by (±£1,2, T^.i)- 
The rectangle (±i 12 , ±^2,1) appears because the chosen 
pairs j( 2 ') do not necessarily have the same I, as was 



the case in (42 1. Interestingly, the (j^ 1 ',^ 2 ') distribution 
is quite non-uniform, but concentrated around (0, 0) and, 
in particular, at the vertices of the rectangles. 

Other mixing angles are difficult to deal with at the 
analytical level, but the (j^, j^) distribution can be ob- 
tained as in the previous case. This is done in Figure [8] 
for ip — it/8, and with the same values of ti and t 2 as be- 



fore. We see that the restrictions (44) and (451 still apply, 
the distribution is depleted around the central point (0, 0), 
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-0.8 -0.4 0.4 0.8 

Fig. 7. Color scale plot of the distribution (j' 1 ', j' 2 ') for t\ = 
0.8 and t% = 0.4 in the case of a left-right symmetric scatterer 
with a channel mixing angle cp — 7r/4. 
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Fig. 8. Color scale plot of the distribution (j' 1 ',! ) f° r ti — 
0.8 and ti = 0.4 in the case of a left-right symmetric scatterer 
with a channel mixing angle p — tt/8. 



and the concentration around the vertices of the rectangles 
is more pronounced than for the case of tp — 7r/4. In Fig- 
ure [9] we present the probability distribution of (j^,]^) 
as a function of the angle (p. Considering different hori- 
zontal cross-sections, like Figures [7] and [8] we can follow 
the evolution from the uncoupled channel case to that of 
intermediate mixings. 



4.2 Asymmetric nanosystems and random channel 
mixing 

In the previous subsection we have seen that in the case 
of left-right inversion symmetry the transmission ampli- 
tudes can be inferred from the study of the distribution of 
matrices M with fixed radial parameters Ai^ and varying 
angular parameters. This procedure can be formalized by 
taking a random matrix approach for the u and v matrices, 
while keeping T fixed. This kind of approach is widely used 
in the application of random matrix theory to quantum 
transport, where we separate the statistics of the matrices 
u and v from that of the transmission eigenvalues. While 
the first ones are assumed to be uniformly distributed 
over the unitary ensemble, the distribution of the latter 
depends on the nature of the problem under study, i.e. 
its diffusive or chaotic character [21]. In practice, such a 
method would require to connect the nanosystem under 
study to mode-mixing reflectionless diffusors, represented 
by matrices uniformly distributed on the unitary ensem- 
ble, and to study the resulting distribution of persistent 
currents. 

Since our procedure will be carried out numerically, we 
can drop the symmetry requirement of the previous sub- 




Fig. 9. Probability density to find a pair (j (1 \ j (2) ) for ti = 0.8 
and £2 — 0.4 with left-right symmetry as a function of the 
mixing angle ip (vertical axis). 



section. We then generate a large number of independent 
unitary matrices u and v for fixed transmission amplitudes 
t\ and £2- The two non-degenerate k- values arising from 
the quantization condition (251 for <P = 7r/2 yield pairs 
of normalized persiste nt c urrents j*- 1 ^ and j^ 2 \ that we 
obtain from Equation (151. The resulting distribution is 
presented in Figure 10 Comparing with Figure [9J there 
are some differences since we do not have spatial symme- 
try and the angle tp is averaged over the unitary ensemble. 
However, similar information can be extracted from both 
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Fig. 10. Color scale plot of the probability density to find a 
pair ft 2 ') for t\ = 0.8 and £2 = 0.4 with random channel 
mixing. 



— ti = 0.8, t 2 = 0.2 
--- t 1 = 0.8, t 2 = 0.4 
tj = 0.8, t 2 = 0.6 
lU " ' ! ii = 0.8,i 2 = 0. 



0.01 



0.001 



Fig. 11. Probability distribution p(j) of the normalized total 
persistent current j in a ring containing a two-channel scatterer 
with random channel mixing for ti = 0.8 and £2 = 0.2 (solid 
line), 0.4 (dashed), 0.6 (dotted) and 0.8 (dashed-dotted) . 



cases: the transmission amplitudes appear as the bounds 
and the most likely values of the j^) distribution. 

If we now consider a non-interacting many-particle sys- 
tem, the sum over the occupied levels should be carried 
out in order to obtain the total persistent current. Assum- 
ing that, like in the N — 1 case, the total persistent cur- 
rent is dominated by the contribution at the Fermi level, 
it is of interest to consider the distribution of the sum 



11 



j = f^+f 2 " 1 ■ Such a distribution is presented in Figure 
on a logarithmic scale for t\ = 0.8 and various values of 



t 2 < t%. The structure found in Figure 10 translates into 
bounds of j at ±(ti + t 2 ) and strong peaks at ±|ii — t 2 \. 
Therefore, from the distribution of j when branching the 
nanosystem to random reflectionless diffusors we can infer 
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Fig. 12. For a two-channel scattering problem there exist six 
different one-channel scattering problems which can be treated 
with the standard embedding method. 



the values of the amplitudes t\ and t 2l and eventually the 
conductance of the system. 

The relationships that we have found between per- 
sistent currents and transmission amplitudes in the non- 
interacting two-channel case are interesting, but difficult 
to use in a practical way. The extension beyond TV = 2 is 
in principle possible. However, as N grows it becomes in- 
creasingly difficult to identify the N transmission eigenval- 
ues from the distribution of j = X^bLi ff b ^ ■ O ur mam inter- 
est is in the N = 2 case once the interactions are switched 
on. The use of the statistical approach for the interacting 
case, leaving practical considerations aside, would rely on 
three key assumptions. Firstly, like in the N = 1 case, we 
have to accept that the interacting problem is described 
by effective one-particle scattering parameters. Secondly, 
the total persistent current should be dominated by its 
Fermi level contribution. Lastly, the two channels at the 
Fermi level need to be quasi-dcgcncratc. 



5 Reduction of multi-channel scattering to 
single-channel problems 

As we have seen in the previous section, it is not obvious 
how to extract the transmission probabilities of an N- 
channel scatterer required for the calculation of the con- 
ductance from the persistent current distribution of an 
Af-channel ring. However, it is possible to calculate sev- 
eral single-channel persistent currents which then allow 
in principle to characterize completely the many-channel 
scattering matrix. 

The effective one-channel scattering problem ij is ob- 
tained by connecting two channels i and j which both 
can be situated on either side of the scattering region. All 
other channels are closed by reinjecting the outgoing cur- 
rent into the scattering region. For the example of two 
channels we present in Figure [12] the six different possible 
setups. Note that the noninteracting ring does not nec- 
essarily connect opposite sides of the interacting region. 
In principle, the electrons can acquire an arbitrary phase 
upon injection. While this phase may represent a useful 
tool, we will set it to zero in the sequel. For an Af-channel 
problem one can construct in this way 2A^ 2 — A^ different 
one-channel scattering problems. 

Within the embedding method, the effective one- 
channel setup linking channels i and j has an associ- 
ated persistent current from which the modulus of the 
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corresponding transmission amplitude can be determined. 
Clearly, the are functions of the 2N 2 + N parameters 
defining the matrices u, v, and A. Such functions are pre- 
sented in the appendix for the two-channel case. In the 
general case, we have 2N 2 — N different T^, so that 2N 
irrelevant phases are left undetermined. 

Our goal is to extract the radial parameters A a that 
determine the conductance. Already the application of the 
embedding method to all of the 2N 2 —N one-channel prob- 
lems may require considerable numerical resources and the 
procedure might be quite lengthy. In addition, the com- 
plicated functional dependence of the makes obtaining 
the A a a challenging task. Depending on the symmetries 
of the system under investigation, both the number of pa- 
rameters appearing in the transfer matrix and the number 
of noncquivalent one-channel problems may be reduced 
with respect to these general expressions. 

Furthermore, interaction-induced nonlocal effects im- 
pose some restrictions. For the effective transmission of 
interacting one-dimensional scatterers, it has been found 
that the usual composition law of scatterers in series does 
not hold [15, 17]. The deviations can be interpreted as a 
nonlocal interaction-induced effect on the effective trans- 
mission of an interacting scatterer. This effect disappears 
only in the limit of large distance between the scatter- 
ers. In the present situation, the closing of channels might 
represent a perturbation leading to such a nonlocal effect 
unless the connecting leads are sufficiently long. 

Despite these difficulties the proposed method might 
provide a possible approach for the demanding problem of 
many-body transport with mixing channels. For instance, 
we might study larger Hubbard chains than the ones we 
have addressed with the other methods. In addition, the 
fact of separating the problem into single-channel setups 
can be useful for the general analysis of the multi-channel 
scattering case as will be exemplified in Appendix [B] 

6 Conclusions 

The embedding method for extracting the conductance of 
an interacting nanosystem is by now a well established 
and widely used procedure in the case of one-dimensional 
spinlcss fermions [9-15]. Its validity has been ascertained 
and interesting physical phenomena have been found from 
its application [14,15,18]. The extension of the embedding 
method to the case of two channels is therefore of great 
importance. In this work we have presented various possi- 
ble extensions of the N = 1 embedding method to N — 2, 
depending on the nature of the channels. 

When the channels represent the two spin directions, 
we have a one-dimensional Hubbard chain, where the ab- 
sence of an effective spin-flip term allows for a straight- 
forward generalization of the N — 1 embedding method. 
The use of the DMRG algorithm in the presence of spin 
makes the problem numerically demanding, but we have 
checked that the method is reliable for chains of up to 4 
interacting sites. We have introduced an improvement of 
the standard embedding method by considering damped 
boundary conditions, allowing to address larger systems 



and to deal with resonances, which is the case for an odd 
number of interacting sites Lg. Our results are compati- 
ble with the perfect conductance predicted for a Hubbard 
chain with an odd number of sites [28]. For even L$ the 
conductance is reduced by the electron-electron interac- 
tions. Similar qualitative behavior has been established 
for Hubbard chains coupled to reservoirs using NRG [29] . 

The generalization of the embedding method to the 
case of two mixed channels is difficult due to the in- 
creased number of transmission amplitudes that are nec- 
essary to determine the conductance. In the case of quasi- 
degenerate channels, by connecting different pairs of chan- 
nels we reduce the problem to various one-dimensional 
setups, where the corresponding transmission amplitude 
can be obtained from the associated persistent current. 
The implementation of such a method requires consider- 
able numerical resources and solving a system of eight 
coupled non- linear equations. In the non-interacting case 
with quasi-dcgcncratc channels we proposed a statistical 
approach that leads to the values of the transmission am- 
plitudes from the distribution of persistent currents ob- 
tained while varying the mixing of the channels. The prac- 
tical implementation of this method in order to extract the 
conductance of a given nanosystem would rely in coupling 
it with an ensemble of reflectionless diffusors that scram- 
ble the phases of the transfer matrix without affecting the 
transmission amplitudes. 

The embedding method is based on an effective one- 
particle description of a correlated system. Anderson-like 
quantum impurity models like that of a short Hubbard 
chain scattering conduction electrons or the Interacting 
Resonant Level Model [4-6] have been studied in the 
literature through the numerical renormalization group 
(NRG) [3,29] or the Bethe-Ansatz [4-6]. From these stud- 
ies we know that the low-energy excitations behave like 
free fermions in the zero-temperature limit of those mod- 
els. This is the case in regimes where the effect of in- 
teractions remains perturbative, as well as in the non- 
perturbative limit. In this latter limit, where (spin or or- 
bital) Kondo effects occur, the extended conduction elec- 
trons continue to behave as free fermions below the Kondo 
temperature, as first shown by Wilson [3]. 

Extending this concept to arbitrary interacting 
nanosystems, one can justify the assumption that an in- 
teracting region embedded between non-interacting leads 
can still be described at sufficiently low temperatures 
and sufficiently far away as an effective one-body scat- 
terer with interaction-dependent effective parameters. In 
certain studies of quantum impurity models, the effec- 
tive quantum transmission is extracted directly from the 
knowledge of the effective one-body excitations or, via the 
Friedel sum rule, from the knowledge [32, 33] of the occu- 
pation numbers of the interacting region. 

Although, not surprisingly, many-channel problems re- 
quire a larger numerical effort, the extension of the em- 
bedding method to more than one channel is feasible and 
we have presented several ideas to achieve this generaliza- 
tion. The embedding method is thus available for future 
practical applications beyond the one-channel case. 



12 



A. Freyn el al: The embedding method beyond the single-channel case 



We thank R.A. Molina for useful discussions and a careful 
reading of the manuscript and H.U. Baranger for clarifications 
about the symmetries in the Hubbard problem. R.A.J, thanks 
the Institute for Nuclear Theory at the University of Washing- 
ton for its hospitality and the Department of Energy for partial 
support during the completion of this work. Financial support 
from the European Union within the MCRTN program and 
from the ANR (grant ANR-08-BLAN-0030-02) are gratefully 
acknowledged. 



A Formulas for the two-channel problem 

In this appendix we present the explicit form of Tr{M} 
and ^(M), which are necessary for solving the quantiza- 
tion condition (23 1 of the two-channel problem. We are 



interested in expressing them in terms of k\ ± k%. For the 
former we have 



Tr{M} = T 1/2 [cos(fc 1J L L - cos(i?) 

+ cos(/c 2 £l - £2) sin(J7)] , 



with 



T 
4 



cos(2</j) cos(2ip) 
2 

1 



(46) 



(47) 



tit 



1'2 



sin(2<,5) sin(2V>) cos(/3 + 77) 



tan£i 



ti tan(<^) tan(^) cos(a + (3 + e + 77) — t 2 cos(a + e) 
t\ tan((/?) tan( , 0) sin(a + (3 + e + 77) — t% sin(a + e) 

(48) 



tan £2 = 

t\ cos(7 + (3 + 9 + 77) — t 2 tan(yj) tan(^) cos(7 + 9) 
ti sin(7 + (3 + 9 + 77) — t 2 tan(</?) tan(^) sin(7 + 9) 

(49) 

tan 2 fl = 

t\ + t\ tan 2 (c/?) tan 2 (-0) — 2t\t 2 tan(^) tan(^) cos(/3 + 77) 
t\ tan 2 (<^>) tan 2 (?/>) + t\ — 2t\t 2 tan(tp) tan(^) cos(/3 + 7/) 

(50) 



The evaluation of (24 1 yields 



F(M) 



+ 



1 + cos(2ip) cos(2-0) 
2 

sin(2</?) sin(2f/>)[cos(/? + 77) — r\r 2 cos(/3 — 77)] 
2~t^t 2 

cos([fci + k 2 ]L L + a + /3 + 7 + e + ?/ + 6») 

hh 

sin(2(^) sin(2^>) cos([fci — k 2 ]Li J + a + e — 7 — 



- ( cos([fci - k 2 ]L L + a + e - 7 



2tit 2 

x [(1 + cos(2</?) cos(2?/>)) cos(/3 + 77) 
+ rir 2 (l — cos(2y>) cos(2i/0) cos(/3 — 77)] 
+ sin([fci — k 2 ]L L + a + e — 7 — 9) 

x [(cos(2^) + cos(2-0)) sin(/3 + 7/) 
+ rir 2 (cos(2ip) — cos(2t/>)) sin(/3 — 77)]) . 




Fig. 13. A Hubbard chain is closed to a ring by a single- 
channel lead attached to different spin channels. Within the 
embedding method, one therefore has to introduce a spin-flip 
in the noninteracting lead. 



The above expressions simplify considerably in the 
case of left-right inversion symmetry, allowing to under- 
take the analytic calculations of Section [4] However, this 
approximation leaves aside important cases, like the one 
of disordered systems. 



B Alternative proof of the spin conservation 
using the reduction of multichannel 
scattering to single-channel problems 



We now demonstrate that the results (|3T| and (32) can 



alternatively be obtained by means of the single-channel 
scattering problems introduced in Section [5] To this end, 
we show that in four of the six scattering configurations 
the persistent current vanishes. 

The number of electrons N a — ^ { fii^ a having spin a 
commutes with the Hubbard Hamiltonian (27 1. Therefore, 
in a closed system the number of electrons having spin up 
or spin down cannot change. While the total spin of the 
electrons inside the interacting nanosystem can fluctuate 
when electrons move between the leads and the nanosys- 
tem, in a situation of stationary transport through the 
system, the mean spin of the electrons entering the sys- 
tem has to equal the mean spin of the electrons leaving 
the nanosystem. It is then possible to describe the effective 
stationary transmission through the Hubbard chain with- 
out invoking transmission processes that are accompanied 
by a spin flip. 



The absence of spin- flip transmissions (31 ) can be con- 
firmed using the formalism presented in Section [5] Four 
of the six possible single-channel configurations connect 
channels of different spin orientation. The spin has then 
to be flipped somewhere in the noninteracting lead as is 
indicated in Figure [T3| for one of these configurations. The 
contribution of the kinetic energy to the total Hamiltonian 
( 27 1 in this example reads 



H K = 



Ls-l L/2-1 

C »,T C i+l,T + C i,l C i+l, 
U=-£/2 »=1 



+ C -L/2,T C £/2,1 + 



(52) 



(51) 



where the third term on the right-hand side is responsible 
for the spin flip. 
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In a single-channel configuration linking different spin 
channels, no persistent current can occur. A non- vanishing 
persistent current in the presence of a spin flip term would 
result in an accumulation of one spin orientation in the 
ring. In a stationary situation, this is of course impossible. 

It can be seen from the relations presented in Ap- 
pendix [C] between the persistent currents in different 
single-channel rings and the effective transmission ampli- 
tudes that the vanishing of persistent currents in setups 
connecting different spin channels is consistent with the 
van ishing of the transmission amplitudes with spin flip 

pil. 



C Reduction of a two-channel scattering 
problem to effective one-channel scattering 
problems 

As illustrated in Figure [12"| for two channels, the effective 
one-channel scattering problem is obtained by closing on 
themselves one pair of channels and connecting the re- 
maining pair by a non-interacting lead. The closing of a 
channel allows in ^ to make the corresponding in- and 
outgoing amplitudes equal, thus reducing by two the num- 
ber of amplitudes involved in the scattering problem. 

The closing of a channel implies that the amplitudes 
of the corresponding in- and outgoing channels become 
equal. Therefore, the closing of two channels reduces the 
number of amplitudes appearing in Q from 8 to 6. Fur- 
thermore, we can use two of these equations in order to 
eliminate the remaining two amplitudes related with the 
closed channels. This leaves us with two equations relating 
the amplitudes of the open channels and thus defining the 
effective one-channel problem. 

This procedure results in a transfer matrix of consid- 
erable complexity when expressed in terms of the param- 
eters of the full scattering problem. The complexity also 
becomes apparent when one expresses the transfer matrix 
in terms of the transmission and reflection coefficients as 
well as the phases encountered on all internal scattering 
sequences leading from one open channel to the other. 

Two different kinds of one-channel setups must be dis- 
tinguished, depending on whether they connect channels 
on opposite sides or on the same side. For the first group 
the moduli of the transmission amplitudes are 



T12 = Qi(~/,6,0,ri,<p,ip) 
T34 = Qi(a, e, (3,Tj,cp- it/2, tp - tt/2) 
T13 = Qi(j, e, (i, 77, ip,ij}- tt/2) 
T 24 = Qi(a, 9, (3,r),(p- tt/2, tp) , 



(53) 



where 



Ql(wi, 07 2 , 07 3 , O7 4 , 07 5 , 07 6 ) 



I N 1 (u) 1 ,uj 2 ,uj 3 ,uj 4 ,uj 5 ,uj 6 ) 2 

L>l(o7i,07 2 , 073,074,U7 5 , U7 6 ) 



(54) 



while for the second group we have 

Tu = Q2(^,^,e,9,r],P) 
T23 = Q2(ip,(p,a,^,f3,ri) 

with 



(55) 



^j N 2 (uJi,u}2, 073, ^4, w 5 ,Q7 6 ) 2 



D 2 (uJ2, 073,075,074) 

The fact that each of the two groups is characterized by a 
single function, Qi or Q 2 , stems from the symmetries of 
the problem. 

In ( 54 1 , we have introduced the definitions 



A^l(07l,O72,O7 3 ,074,O7 5 ,O7 6 ) 

! 1 COS(075) COs(076)[r2 COs(o7i + 073 — 072 — W4) 

— COS(071 + 07 3 + 07 2 + 07 4 ) 
— t2 Sin(075) Sin(o7g)[r! COs(o7i — 07 2 ) 

— COS(07! + 072)] 



and 



(57) 



Dl{uJl, 07 2 ,O7 3 , 07 4 , 075 , 07 6 ) = 

[7Jl(o7l,07 2 ,073,O74 , 075,O7 6 ) - l] 2 

+ #l(wi,072 , 073,O74 , 07 5 ,O7 6 ) 2 (58) 

where 

^(071,072,073 , 074 , 075 , 076) = 

ri J 4 2 (o7i,o7 2 , 075 - 7r/2,076 - tt/2) 
+r 2 A 2 (uj 1 + 073, o7 2 + ^4, W5, we) 
-rir 2 A 3 (uJi + 07 2 + W3, 07i + 07 2 + w 4 , 075, 07 6 - tt/2) 

+7^(071 + 07 2 , 07 3 + 07 4 , 07 5 - 7r/2 , 07 6 - 7r/2) 
— A\{UJX + 07 2 + W3 + 07 4 , 07l + 07 2 , 07 5 , 07 6 ) 



with i = c, s and 



(59) 



^l(wi,07 2 ,073,07 4 ) = 

2*1*2 COS(073) COS(074) sin(o73) sin(o74) COs(2o7i + 07 2 ) 

^1(071,072,073,074) = 

2*1*2 cos(o73) cos(o74) sin(o73) sin(o74) sin(2o7i + 072) 
A 2 (uji,lu 2 , 073, 074) = cos(o7 3 ) 2 cos(2o7i) 

+ COS(074) 2 COS(2072) 

A s 2 (071 , 07 2 , 073 , 074 ) = cos(o7 3 ) 2 sin(2o7i ) 

+ cos (074) 2 sin (2072) 

Al(ldi,UJ 2 , 07 3 , 07 4 ) = COS(073) 2 COS(074) 2 COs(2o7i) 

+ sin(o7 3 ) 2 sin(o74) 2 cos(2o7 2 ) 
A 3 (lui, uj 2 , 073, 074) = cos(o73) 2 cos(o74) 2 sin(2o7i) 

+ sin(o7 3 ) 2 sin(o74) 2 sin(2o72) . 

(60) 
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The functions appearing in Q 2 are defined as 16. A.O. Gogolin, N.V. Prokof'ev, Phys. Rev. B 50, 4921 

(1994) 

N2(u}i,0J2,u>3,uj4,uj 5 ,uJe) = 17. D. Weinmann, R.A. Jalabert, A. Freyn, G.-L. Ingold, J.-L. 

ht 2 cos(2 Wl ) sin(2 W2 ) sin(^ 3 - w 4 ) Pichard > Eur ' Ph y s ' J ' B 66 ' 239 ( 2008 ) 

18. G. Vasseur, Ph.D. thesis, Universite Louis Pasteur 
-sin(2wi)(sin(w 2 ) sin(w 3 - w 4 + uj 5 + uj 6 ) Strasbourg, 2006, e-pri nt |http: //eprints-scd-ulp~| 



-COs(u; 2 ) 2 sin(w 3 - W 4 - C0 5 - U e ) u-strasbg.fr : 8080/574/ 

(61) 19. PA. Mello, J.-L. Pichard, J. Phys. I 1, 493 (1991) 
20. C.W.J. Beenakker, Rev. Mod. Phys. 69, 731 (1997) 



-r<i sin(w 3 + UJ4 + uj 5 — uu§) 



-ri sin(w 3 + LUi + lu 5 + Lu e ) 21. R.A. Jalabert, J.-L. Pichard, J. Phys. I 5, 287 (1995) 



-r 4 r 2 [cos(w 2 ) 2 sin(w 3 — w 4 — uj 5 + uj 6 ) 



22. X. Waintal, G. Fleury, K. Kazymyrenko, M. Houzet, 
P. Schmitteckert, D. Weinmann, Phys. Rev. Lett. 101, 

- sin(cj 2 ) 2 sin(cj 3 - w 4 + uj 5 - u 6 )]) , 106804 (2008) 

23. S.R. White, Phys. Rev. Lett. 69, 2863 (1992) 

_D 2 (c<j 2 , (jj 3 , UI5, UI4) — [i? 2 (ijJ 2 , CiJ 3 , [JJ5, w 4 ) — l] 2 24. Density-Matrix Renormalization — A New Numerical 

(62) Method in Physics, edited by I. Peschel, X. Wang, M. 



if 2 (w 2 , w 3 ,w 5 ,u 4 ) = riA 4 (w 2 , w 3 ,w 4 ) 



Kaulke, K. Hallberg, Lect. Notes Phys. 528, 1-355 
(Springer, Berlin, 1999) 
25. P. Schmitteckert, Ph.D. thesis, Universitat Augsburg, 1996 
+ r 2 A 4 (w 2 ,w 2 + cj 4 ,lji + u s ) 26. A. Oguri, Phys. Rev. B 63, 115305 (2001) 

- rir 2 cos[2(w 3 + uj 5 + lu 4 )} 27. Y. Nisikawa, A. Oguri, Phys. Rev. B 73, 125108 (2006) 

(63) 28. A. Oguri, Phys. Rev. B 59, 12240 (1999) 

29. A. Oguri, A.C. Hewson, J. Phys. Soc. Jpn. 74, 988 (2005) 
+ r 2 A s 4 (u;2,uj2 + uj 4 ,uJi + uj 5 ) 3 q p. Schmitteckert, R. Werner, Phys. Rev. B 69, 195115 

- r x r 2 sin[2(w 3 + u 5 + lu 4 )} , (2004) 

31. D. Bohr, P. Schmitteckert, P. Woffle, Europhys. Lett. 73, 
and 246 (2006) 



#2(0/2, w 3 ,o;5,w 4 ) = riA\(uj 2 , ^3,w 4 ) 



2„„„/o x t r . ( ^ 2 _ /r> s 32. T.K. N'g, P.A. Lee, Phys. Rev. Lett. 61, 1768 (1988) 

33. A. Freyn, J.-L. Pichard, e-print |arXiv:0904 .0889^1 



A1(uj2, ^3, w 4 ) = cos(w 2 ) cos(2lj 3 ) + sin(w 2 ) cos(2o; 4 ) 

A%(u>2, ^3, — cos(w 2 ) 2 sin(2w 3 ) + sin(w 2 ) 2 sm(2w 4 ) . 

(64) 
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